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APPLICATION  OF  ADAPTIVE  CONTROL 
TO  AIRCRAFT  SYSTEMS 

Chen  Zongji  and  Zhang  Hongyue 
( Beijing  Institute  of  Aeronautics  and  Astronautics ) 


Abstract 

^his  paper  reviews  the  competition  of  the  adaptive  technique  with 
other  techniques.  The  inherent  drawback  of  the  adaptive  technique  in 
application  to  aircraft  systems  is  considered  in  the  light  of  the  analysis 
of  the  features  of  flight  control  systems.  The  potentials  and  perspectives 
of  the  adaptive  technique  are  also  suggested. 

The  development  of  advanced  aircraft  demands  enhanced  performances 
which  go  beyond  the  capacity  of  airdata  scheduling  technique.  Therefore, 
the  adaptive  control  schemes  are  expected  to  be  applied  to  these  flight 
control  problems.  In  order  to  make  full  use  of  the  advantages,  such  as 
saving  sensor  hardware,  requiring  less  priori  knowledge  of  aircraft  and 
providing  higher  flexibility,  it  is  necessary  to  carry  out  robustness  re¬ 
search  of  existing  adaptive  control  schemes  and  to  develop  new  more  ro¬ 
bust  adaptive  control  schemes.*N  v  \J\  J 
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In  1979,  Stein^1^,  Kreisselmeier  ^ 2  ^  and 
.  r  3  I 

Rynaski  presented  separate  papers  on  the 
application  of  adaptive  control  to  aircraft 
systems  in  the  symposium  for  the  Application 
of  Adaptive  Control.  This  paper  will,  through 
a  general  review  of  these  papers,  establish 
a  clearer  knowledge  on  the  application  status 
of  adaptive  aircraft  control,  the  problems  it 
faces  and  its  possible  application  direction 
and  perspectives. 


I.  Preface 


For  twenty  or  more  years,  many  theoreticians  and  practical  engi¬ 
neering  technicians  have  expended  much  effort  trying  to  apply  the 
adaptive  control  to  aircraft  systems.  However,  very  few  results 


have  been  obtained  up  to  now  despite  the  efforts. 


Theoretically,  the  aircraft  control  problems  which  have  the 
characteristics  of  parameter  uncertainty  and  slow  time-variance  should 
be  a  perfect  domain  for  the  application  of  the  adaptive  control. 

Yet,  why  has  the  application  of  adaptive  control  to  aircraft  systems 
been  so  limited?  This  is  a  question  that  has  puzzled  many  adaptive 
control  theoreticians  and  aeronautical  engineering  technicians  who 
are  interested  in  adaptive  control.  In  1979,  in  the  Symposium  for 
the  Application  of  Adaptive  Control  held  by  Yale  University,  Stein, 
Kreisselmeier  and  Rynaski  presented  papers  on  the  application  of 
adaptive  control  to  aircraft  systems.  These  were  more  encompassing, 
more  in-depth  papers  with  substances  prepared  by  people  who  tried 
to  answer  the  above  question  in  recent  years.  This  paper  will,  on 
the  basis  of  introducing  these  papers  and  combining  with  the  develop¬ 
ment  of  adaptive  control  in  recent  years,  conduct  a  general  review 
on  the  application  status  of  adaptive  control  to  aircraft  systems, 
existing  special  problems,  its  possible  application  direction  and 
perspectives. 

II.  Features  of  an  Aircraft  and  Its  Control  Systems 

The  rigid  aircraft  is  a  moving  object  of  the  complex  six  degrees 
of  freedom  ana  multiple  maneuvering  surfaces.  Its  demands  for  high 
speed,  high  mobility  and  high  reliability  make  the  control  problem 
of  the  aircraft  more  complicated.  In  addition,  the  aircraft  control 
systems  have  the  following  two  special  problems,  i.e.,  the  uncertainty 
of  the  controlling  object's  parameter  and  model,  and  the  man-aircraft 


1.  Uncertainty 


The  uncertainties  of  aircraft  are  expressed  in  two  areas,  i.e., 
parameter  uncertainty  and  model  uncertainty. 

The  parameter  uncertainties  of  aircraft  are  caused  by  the  flight 

conditions  with  a  large  range  of  change  and  the  change  in  aircraft 
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structure.  The  change  in  flight  conditions  mainly  means  the  change 
in  flight  altitude,  flight  velocity  and  flight  posture.  The  change 
in  aircraft  structure  mainly  means  the  change  in  fuel,  loads  and 
geometric  dimensions  during  the  flight  process.  These  changes  cause 
the  parameters  in  the  aircraft's  dynamics  and  kinematics  equations 
to  change.  During  the  flight  process,  the  changes  are  actually  un¬ 
predictable.  However,  research  has  shown  that  numerous  parameter 
changes  can  be  grouped  eventually  into  a  few  parameter  changes  of 
slow  time-variance.  Take  the  vertical  control  problem,  for  example: 
the  dominant  parameters  are  the  pitching  moment  Mg  generated  by 

the  stabilizer  and  the  pitching  moment  M#  generated  by  the  attack 
[  4  1 

angle  .  It  must  be  pointed  out  that  the  flight  enveloping  lines 
of  modern  aircraft  have  a  tendency  to  increase  to  a  large  attack 
angle  and  post-stall  ability.  However,  we  still  lack  sufficient 
knowledge  on  the  aircraft  dynamic  responses  under  this  kind  of  flight 
condition.  Moreover,  the  change  of  these  parameters  are  no  longer 
of  slow  time-variance. 

The  uncertainty  of  the  model  means  the  unknown  high  frequency 
dynamic  responses  or  those  that  are  not  considered  during  the  process 
of  model  building.  They  include  high  frequency  structural  harmonic 


vibration,  unstable  air  distribution  dynamic  responses,  hydraulic 


mechanism,  transfer  delay  of  aerodynamic  responses,  etc.  In  essence, 


aircraft  is  a  dynamic  system  of  infinite  dimension,  but  low  order 


systems  are  generally  used  to  approximate  in  designing.  For  example, 


the  low  frequency  input-output  responses  of  the  aircraft  in  vertical 


motion  can  be  approximated  by  a  third-order  rigid  model  or  an  elastic 


model  with  a  little  higher  order,  but  there  is  great  difference  between 


its  high  frequency  end  and  the  third-order  model.  Figure  1  gives 


the  Bode  diagram  of  the  third-order  transfer  function.  Meanwhile, 


it  also  gives  the  enveloping  lines  of  the  actual  aircraft  frequency 


response.  The  frequency  response  of  the  aircraft  high  frequency 


end  can  be  any  continuous,  smooth  curve  in  the  shaded  area  of  the 


f  igure . 


( 2 ) 


Fig.  I . 


frequency  response  and 


Key:  (lj  Gam;  (2)  Frequency;  (3)  Phase. 


requency  uncertainties 


2.  Man-Aircraft  System  and  Flight  Quality 


The  purpose  of  any  flight  control  system  is  to  improve  or  increase 


the  flight  quality  of  aircraft,  and  the  flight  quality  defines  the 


coordinating  relations  of  the  man-aircraft  system.  The  flight  quality 


of  the  vertical  short  cycle  motion  is  related  to  the  three  parameters 
of  the  aircraft  dynamic  response,  and  they  are  the  zero  point,  inherent 
frequency  and  damping  ratio  of  the  pitching  speed  with  respect  to 
the  stabilizer  transfer  function. 

According  to  the  U.S.  military  aircraft  flight  quality  specifica- 
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tion  MIL-F-8785B( ASG) L  ,  for  first  class  flight  quality  standard', 
the  allowable  change  for  the  damping  ratio  is  about  4  times  and  the 
allowable  change  for  the  inherent  frequency  is  about  10  times.  Accord¬ 
ingly,  the  pilot  has  strong  adaptability  to  the  large  range  change 
in  aircraft  inherent  frequency,  i.e.,  the  man-aircraft  system  actually 
includes  an  "adaptive  controller " -pi lot . 

Yet,  research  has  shown  that  the  pilot's  adaptability  to  time 
delay  is  weaken.  When  the  pilot  gives  an  instruction  to  the  aircraft, 
he  expects  the  aircraft  to  have  immediate  response  to  the  instruction. 
If  there  is  a  100  ms  delay  caused  by  the  computation  time  of  the  com¬ 
puter  and  other  factors,  this  will  seriously  affet  the  pilot.  In 
addition,  the  pilot's  adaptability  to  the  response  of  the  nonminimal 
phase  is  alio  very  poor.  The  nonminimal  phase  will  cause  its  initial 
response  to  be  in  the  opposite  direction  of  the  demand  of  the  instruc¬ 
tion. 

From  the  standpoint  of  flight  quality  demand  of  the  man-aircraft 
system,  the  dynamic  response  of  the  aircraft  is  allowed  to  vary  within 
a  larger  range,  but  the  higher  order  dynamic  delay,  transfer  delay 
of  the  numerical  computer  and  the  nonminimal  phase  should  be  avoided 


III.  Adaptive  Control  and  Its  Control  Scheme  in  the  Competition 
History  of  Aircraft  Systems 

In  the  history  of  aeronautical  development,  the  competition  among 
various  control  schemes  has  usually  been  fierce.  People  generally 
think  that  the  adaptive  technique  performs  better  than  the  regular 
feedback  system  in  systems  with  uncertain  situations.  Therefore, 
from  the  stand  point  of  uncertainty  of  the-aircraft  dynamic  response, 
the  adaptive  technique,  in  principle,  should  be  a  more  ideal  control 
scheme.  However,  whether  the  aircraft  manufacturing  companies  and 
customers  are  willing  to  adopt  this  scheme  still  requires  considera¬ 
tion  of  aircraft  performance,  cost  effectiveness,  reliability,  etc. 
provided  by  said  scheme.  In  the  early  60 ' s  and  mid  70's,  there  had 
been  two  rounds  of  competition  between  the  adaptive  technique  and 
other  control  schemes,  and  each  time  was  ended  with  the  defeat  of 
the  adaptive  control  scheme. 

Lne  primary  competitor  of  the  adaptive  control  scheme  is  the 
airdata  scheduling  control  laws.  It  can  be  seen  from  section  II  that 
the  change  of  aircraft  parameters  can  be  grouped  into  the  change  of 
certain  dominant  parameters,  such  as  and  ,  and  these  dominant 

parameters  are  closely  related  to  airdata  kinetic  pressure,  Mach  num¬ 
ber  or  attack  angle.  Thus,  the  control  laws  expressed  in  M^  and 
Ma  can  be  approximated  by  the  control  laws  expressed  in  kinetic 
pressure,  Mach  number  and  attack  angle. 

As  early  as  the  60's,  people  had  higher  requirements  for  flight 
enveloping  lines  and  performance  of  military  aircraft.  But  the  manu- 


facture  and  maintenance  of  an  airdata  sensor  with  high  accuracy  and 
good  reliability  at  the  time  could  not  satisfy  the  requirements  for 
the  airdata  scheduling  control;  thus,  the  adaptive  control  scheme 
was  once  a  very  popular  scheme  and  was  involved  in  the  test  flights 
of  the  X-15,  F-lll,  F-101,  and  F-4  test  aircraft.  The  test  flight 
results  were  not  good,  and  especially  the  crash  of  the  X-15  test  air¬ 
craft  had  a  bad  psychological  influence  on  the  application  of  adaptive 
control  to  aircraft  systems.  At  the  same  time,  with  this,  the  accuracy 
and  maintainability  of  the  airdata  sensor  had  been  improved,  and  its 
cost  had  dropped.  Consequently,  airdata  scheduling  control  laws  became 
the  primary  scheme  for  various  aircraft  control  systems  thereafter. 
Therefore,  some  people  think  that  the  control  engineer  had  lost  to 
the  sensor  manufacturer  in  the  first  round  of  competition  between 
the  adaptive  control  scheme  and  the  airdata  scheduling  control  scheme. 

However,  the  competion  between  the  two  schemes  did  not  end. 

In  the  mid  70's,  another  generation  of  advocates  for  adaptive  control 
started  the  second  round  of  competion,  and  their  contentions  are: 

1.  There  have  been  new  developments  in  the  theory  of  adaptive 
control,  and  they  are  primarily  displayed  in  the  areas  of:  development 
in  the  identification  theory  (e.g.,  parameter  estimation  using  on¬ 
line  maximum  approximation  method);  development  in  the  stability  analy¬ 
sis  and  verification  of  adaptive  control  (the  use  of  [Livapnov’s] 
stability  and  superstabilitv  theories  to  design);  further  perfection 
of  the  theory  and  application  of  the  self-tuning  regulator. 


a 


2.  The  rapid  development  in  computer  software  and  hardware  tech¬ 
nologies  make  some  of  the  complex  adaptive  control  laws  easy  to  be 

implemented.  For  example,  the  F-8  on-board  computer  AP-101  needs 
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only  5  ms  m  every  sampling  period  to  complete  the  solution  of 
on-line  parameter  estimation  and  control  laws  using  the  maximum  approx 
imation  method  when  the  sampling  frequency  is  50  Hz. 

3.  The  adoption  of  a  fly-by-wire  control  sys'tem  must  require 
the  use  of  the  redundancy  technique.  Consequently,  it  is  inevitable 
that  the  circuitry  of  the  airdata  sensor  and  measurement  will  multiply 
manyfold  in  the  scheduling  control  laws.  This  kind  of  increase  not 
only  directly  raises  the  cost,  but  also  causes  difficulties  as  to 
finding  proper  places  to  install  these  sensors. 

4.  some  of  the  control  laws  cannot  be  simply  expressed  by  kinetic 
pressure,  Mach  number  or  attack  angle  (e.g.,  the  active  control  and 
direct  force  control  modes -of  oscillation/tremor). 

The  above  contentions  received  much  attention  from  the  aircraft 
control  industry.  NASA  had  arranged  adaptive  control  testing  in  the 
digital  fly- by -wire  control  system  experimental  project  of  F-8C. 

But  the  testing  results  showed  that  the  performance  of  aircraft  employ 
ing  the  adaptive  control  scheme  was  not  as  good  as  expected;  it  did 
not  even  exceed  the  performance  of  the  airdata  scheduling  control 
scheme.  The  adaptive  scheme  used  in  the  F-8C  and  based  on  the  on¬ 
line  identification  of  Mg  ,  M$<  and  air  velocity  may  perform  satisfac¬ 
torily  only  under  the  effects  of  sustained  and  excited  experimental 
signals.  Therefore,  the  adaptive  control  scheme  represented  by  the 


self  tuning  adaptive  scheme  had  lost  again  in  the  second  round  competi¬ 
tion  with  the  airdata  scheduling  control  scheme. 

IV.  Existing  Problems  in  the  Application  of  the  Adaptive  Technique 
to  Aircraft  Systems 

This  section  will,  through  the  analysis  of  the  existing  problems 
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in  the  application  of  adaptive  control  to  aircraft  systems,  attempt 
to  understand  why  the  adaptive  control  scheme  had  failed  in  past  com¬ 
petitions  . 

The  two  primary  schemes  of  the  adaptive  control  are:  direct 
adaptive  control  and  indirect  adaptive  control.  The  direct  adaptive 
control  is  also  known  as  the  model  reference  adaptive.  This  scheme 
does  not  need  to  conduct  identification  of  the  system's  parameters; 
it  adjusts  the  controller  parameters  according  to  a  certain  stability 
standard,  making  the  output  of  the  controlled  object  to  follow  asymp¬ 
totically  the  output  of  .the  reference  model.  Contrary  to  this,  the 
indirect  adaptive  control,  represented  by  the  self  tuning  regulator, 
first  needs  to  conduct  identification  of  the  system's  parameters, 
then ,  based  on  these  parameters  and  according  to  a  certain  standard, 
the  appropriate  control  laws  are  generalized.  We  will  analyze  below, 
respectively,  what  problems  will  occur  when  the  two  schemes  are  applied 
to  aircraft  systems. 


1.  Model  Reference  Adaptive  Control 


The  structural  principle  of  the  model  reference  adaptive  control 
system  is  as  shown  in  Fig.  2.  Basically,  it  uses  the  feedback  com¬ 
pensator  with  adjustable  parameters  to  eliminate  the  object  zero  point 
and  assigned  model  extreme  points;  the  model  zero  point  is  introduced 
using  the  [feedforward]  compensator,  making  the  combined  system  o.f 
the  controlled  object,  [feedforward]  and  fee'dba'ck'  compensators  and 
the  reference  model  of  equal  value.  In  this  process,  since  the  zero 
point  of  the  object  is  eliminated,  the  model  reference  adaptive  scheme 
cannot  be  applied  to  the  nonminimal  phase  system. 


Fig.  2.  1  reference  adaptive  control  system 

Key:  (i)  [Feedforward]  .compensator;  (2)  Reference  model;  (3)  Controll¬ 

ed  object;  (4)  Feedback  compensator. 

We  also  noticed  that  the  parameters  of  these  compensators  are 
regulated  based  on  the  errors  between  the  object  and  model  output, 
thus,  such  control  laws  cause  serious  nonlinearity.  This  will  affect 
the  pilot's  judgment  and  maneuvering. 

The  verification  of  the  asymptotic  stability  of  the  model  refer¬ 
ence  adaptive  system  is  conducted  under  a  series  of  assumptions.  These 
assumptions  are:  no  external  interference,  the  highest  order  and 


the  relative  order  (the  order  difference  of  the  denominator  and  numer¬ 
ator)  of  the  object  are  known,  minimalphase  system  and  the  sign  of 
gain  of  the  object  is  known.  However,  since  these  conditions  are 
not  all  present  in  our  actual  design  problem,  we  still  lack  clear 
knowledge  on  the  stability  of  the  entire  system  when  the  adaptive 
control  scheme  is  applied  to  the  actual  aircraft  control  systems. 

From  the  standpoint  of  the  multi-purposes  "of'  an  aircraft,  a  singl 
reference  model  is  not  capable  of  handling  it.  For  example,  the  dynam 
ic  effects  of  an  aircraft  during  landing  are  impossible  to  be  as  fast 
as  those  at  high  flying  speed.  Consequently,  the  reference  model 
will  be  adjusted  according  to  flight  duty,  and  this  will  increase 
the  complexity  in  design. 

The  advantages  of  the  model  reference  adaptive  are  that  it  does 
not  require  sustained  and  excited  signals,  nor  does  it  require  that 
the  adjustable  parameters  of  the  controller  converge  uniformly,  and 
the  object  output  is  guaranteed  to  follow  asymptotically  the  model 
output . 

It  can  be  seen  from  the  above  analyses  that  before  the  robust¬ 
ness  of  the  model  reference  adaptive  system  for  external  interference 
and  model  error  is  verified,  the  aircraft  manufacturer  and  the  custom¬ 
ers  cannot  accept  this  scheme  psychologically.  Meanwhile,  the  intro¬ 
duction  of  serious  nonlinearity  and  the  nonapl 1 icabi 1 i ty  to  the  non- 
minimal  phase  system  also  bring  about  certain  difficulties  in  the 
application  of  the  said  scheme. 


2.  Self  Tuning  Regulator 


The  structural  principle  of  the  self  tuning  regulator  is  shown 
in' Fig.  3.  It  conducts  on-line  identification  of  the  system  parameters 
first,  then  the  obtained  parameter  estimation  is  treated  as  the  true 
value  of  the  system  parameter  and  the  control  laws  are  combined  accord¬ 
ing  to  different  design  standards.  These  design  standards  can  be. 
model  tracing,  extreme  points  assignment,  minimum'  direction  difference, 
etc.  Therefore,  the  self  tuning  regulator  is  more  flexible. 


Fig.  3.  Self  :^ning  regulator 

Key:  (1)  Controller;  (2)  Design  calculator?  (3)  Identifier; 

(4)  Self  tuning  controller;  (5)  Controlled  object. 

Kreisselmeier  completed  the  verification^7-^  of  full  stability 
under  sustained  and  excited  conditions  for  one  kind  of  adaptive  ob¬ 
server  and  condition  feedback  control  laws.  Yet,  since  the  regulation 
law  of  the  current  self  tuning  regulator  is  not  built  upon  a  stability 
theory,  the  verification  of  general  full  stability  has  yet  to  be  com¬ 
pleted  . 

It  should  also  be  pointed  out  that  the  on-line  parameter  estima¬ 
tion  up  to  now  is  still  not  a  mature  technique.  Even  the  off-line 
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parameter  estimation,  according  to  experience,  is  a  method  whose  accu¬ 
racy  is  difficult  to  guarantee.  It  is  determined  by  the  dynamic  fea¬ 
tures  of  the  object  accuracy  of  the  sensor,  motion  features  of  the 
aircraft  and  the  identification  method  adopted.  Since  the  success 
or  failure  of  this  kind  of  scheme  is  determined  by  the  accuracy  of 
the  parameter  estimation,  this  requires  that  the  object  be  subject 
to  sustained  and  excited  input  signals  to  assure  that  the  parameter 
estimation  value  converges  at  the  true  parameter  value.  Kreisselmeier 
believes  that,  in  the  flight  evnironment  of  interference  caused  by 
scattered  wind,  the  additional  excited  signals  are  no  longer  needed. 
Stein,  however,  points  out  that  the  test  flight  of  F-8C  showed  that 
the  excited  signals  were  the  key  to  obtaining  satisfactory  performance. 

The  only  information  required  before  designing  the  self  tuning 
regulator  is  the  order  of  the  system,  unlike  the  model  reference 
adaptive  which  requires  more  prior  knowledge.  Also,  the  self  tuning 
regulator  designed  according  to  extreme  point  assignment  and  general¬ 
ized  standards  for  minimum  direction  difference  can  be  applied  to 
the  nonminimal  phase  system. 

According  to  the  above  analyses,  it  is  believed  that  the  self 
tuning  adaptive  control  scheme  is  more  robust  is  the  apllication  to 
aircraft  control  systems.  But  the  self  tuning  adaptive  still  lacks 
clear  conclusions  in  the  area  of  full  stability  and  requires  sustained 
and  excited  signals.  This  also  hampers  the  smooth  application  of 


this  scheme. 


V.  The  Possible  Application  Direction  and  Perspectives  of  the  Adaptive 
Control  in  the  Aeronautical  Industry 

Since  the  adaptive  control  cannot  demonstrate  its  superiority 
in  situations  where  the  airdata  scheduling  control  laws  also  apply, 
we  should  search  for  the  applications  of  the  adaptive  contorl  is 
situations  where  the  scheduling  control  laws  do  not  work. 

For  example,  in  the  high  speed  flight  during  re-entry,  the  airdata 
cannot  be  timely  and  accurately  measured;  for  the  control  of  a  small 
guided  missile,  it  is  not  only  ineconomical  to  install  a  highly  sensi¬ 
tive  sensor  and  corresponding  measurement  circuitry  for  scheduling, 
but  also  will  be  limited  by  its  structure;  in  the  active  control 
duties,  such  as  direct  force  control  and  oscillation/tremor, 
etc.,  the  interrelation  between  the  airdata  and  the  parameters  of 
the  controlled  object  is  very  poor.  The  scheduling  control  laws  are 
not  suitable  for  these  situations.  We  will  introduce  in  more  detail 
below  why  the  airdata  scheduling  control  does  not  work  in  direct  force 
control  and  oscillation/tremor  suppression. 

1.  Direct  Force  Control 


The  primary  trend  of  improvement  of  modern  fighter  performance 
is  to  increase  the  so-called  performance  of  direct  force  cotnrol. 

The  multiple  control  surfaces  of  this  kind  of  aircraft  can  use  the 
coordinated  maneuvers  of  a  certain  coupling  form  to  generate  direct 
lift  or  lateral  force,  but  not  the  pitching  or  deviation  moment. 

This  has  greatly  increased  the  performance  of  aircraft.  In  the  desigr 


of  this  kind  of  control  system,  we  must  know  the  aerodynamic  response 
parameters  of  each  control  surface.  Yet,  the  interrelation  between 
these  response  parameters  and  the  airdata  is  rather  poor,  making  it 
impossible  to  adjust  the  controller  parameters  according  to  airdata. 

2.  Oscillation/Tremor  Suppression 

For  light  structure  or  high  speed  aircraft,  the  interplay  between 
the  structural  dynamic  response  and  air  flow  becomes  more  distinctive. 
When  a  light  damping  vibration  is  generated  by  these  interplays,  it 
is  called  oscillation;  when  these  vibrations  are  unstable,  they  are 
called  tremor.  These  phenomena  all  occur  during  high  speed  flight. 

In  order  to  prevent  aircraft  structural  fatigue  and  crash,  the  speed 
of  aircraft  is  somewhat  limited.  In  recent  years,  the  expansion  of 
such  speed  limit  has  gained  extensive  attention.  However,  due  to 
the  complexity  of  the  oscillation/tremor  model,  the  research  in  this 
kind  of  control  problem  is  progressing  quite  slowly.  The  conlexitv 
of  the  oscillation/tremor  model  is  displayed  in  thut  it  includes  numer¬ 
ous  structural  vibration  modes  that  interplay  with  air  flow  in  the 
10^  200  Hz  high  frequency  end,  and  the  uncertainty  of  these  modes 
as  the  flight  conditions  change.  The  active  control  duty  of  oscilla¬ 
tion/tremor  is  to  measure  the  vibration  mode  parameters  that  require 
suppression  at  appropriate  locations  using  the  appropriate  sensor, 
and  then  increase  the  damping  of  these  modes  by  designing  an  appropri¬ 
ate  compensator.  When  it  is  necessary  to  design  the  osc i 1 la t ion/t remor 
controller  in  full  flight  enveloping  lines,  the  location  of  the  sensor 
and  the  parameter  of  the  compensator  must  all  be  adjusted  as  the  mode 
changes.  Since  the  airdata  are  only  partially  related  to  these  modes, 
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they  still  do  not  possess  this  kind  of  adjustment  means  up  to  now. 

The  adaptive  control  will  be  one  of  the  methods  that  could  possible 
solve  this  type  of  problem. 

Stein,  Kreisselmeier  and  Rynaski  also  believe  that: 

(1)  The  self  tuning  adaptive  scheme  is  more  robust  in  the  applica¬ 
tion  to  aircraft  systems  than  the  model  reference  adaptive,  since 

the  model  building  error  of  the  model  reference  adaptive  will  introduce 
higher  order  dynamic  influences. 

(2)  Unless  the  speed  of  the  numerical  computer  is  very  high  to 
keep  the  computation  delay  less  than  50  ms,  it  is  still  better  to 
use  a  simulation  controller. 

(3)  Measures  must  be  taken  to  increase  the  accuracy  of  the  on¬ 
line  parameter  identification,  since  identification  error  can  also 
introduce  high  order  dynamic  delay. 

It  must  be  pointed  out  here  that  their  conclusions  were  based 
on  the  results  in  the  theory  and  application  of  adaptive  control  before 
1979. 


In  recent  years,  the  theoreticians  of  adaptive  control  have  start¬ 
ed  to  pay  attention  to  the  research  of  robustness  of  the  model  refer- 

r  g  I 

ence  and  model  error.  Kreisselmeier  and  Narendra1  1  verified  that 

under  bounded  disturbances  all  the  signals  of  the  model  reference 

r  g  1 

adaptive  system  are  bounded.  Kokotovic  analyzed  that  when  the 


high  speed  mode  has  weak  influences  on  the  system,  or  the  high  speed 
mode  with  strong  influences  is  subject  to  low  general  filtering,  the 
model  reference  adaptive  system  has  robustness  with  respect  to  model 
error.  The  authors  of  this  paper^^'  1 1  ^  analyzed,  using  the  stability 
theory,  that  the  model  reference  adaptive  system  possesses  robustness 
with  respect  to  the  external  interference  and  model  building  error. 
There  has  been  progress  in  how  to  design  a  model  reference  adaptive 

-  •  r  1 2 1 

system  with  strong  robustness.  Peterson  and  Narendra  have  estab¬ 

lished  a  nonlinear  adaptive  controller  to  suppress  the  influences 
on  system  by  external  interference.  At  the  same  time  with  this 
development,  the  verification  of  the  full  stability  of  the  self  tuning 
adaptive  has  not  obtained  a  true  breakthrough,  and  it  seems  that  the 
requirement  for  sustained  and  excited  signals  is  inevitable.  From 
the  above  analyses,  the  point  that  the  model  reference  adaptive  will 
be  more  robust  in  applications  from  now  on  is  not  unfounded. 

Recently,  the  theoretical  study  of  constant  parameter  robust 
control  and  the  research  in  its  application  to  aircraft  control  systems 
have  gained  substantial  progress^^'  ^  ^  .  In  order  to  develop  purpose¬ 
fully  the  study  of  applications  of  various  techniques,  we  must  clearly 
understand  under  what  conditions  can  the  constant  parameter  robust 
controller  satisfy  the  performance  requirements;  under  what  conditions 
will  the  airdata  scheduling  controller  be  a  better  control  scheme ; 
and  under  what  conditions  is  the  adaptive  controller  a  must.  The 
research  in  this  area  requires  concerted  efforts  by  the  theoreticians 
and  the  flight  control  engineers. 


■  hviiium! 


As  of  now,  the  existing  problems  with  the  airdata  scheduling 
control  laws  not  only  still  remain  unsolved,  but  will  also  become 
more  noticeable  as  the  flight  enveloping  lines  of  modern  aircraft 
further  enlarge  and  its  mobility  further  increases.  At  the  same  time 
with  this,  the  further  perfection  of  adaptive  control  theory  and  the 
rapid  development  of  computer  software  and  hardware  technologies 
will  make  the  application  of  adaptive  control  to  flight  control  systems 
more  and  more  mature.  All  these  indicate  that,  in  the  development 
of  aircraft  control  systems,  the  third  round  competition  between  the 
adaptive  control  and  the  airdata  scheduling  control  schemes  could 
be  inevitable. 

VI .  Conclusions 


Due  to  the  speciil  features  and  inherent  defects  of  the  aircraft 
control  systems,  the  application  of  adaptive  control  to  flight  control 
systems  is  not  successful  up  to  now.  However,  as  the  performance 
requirements  of  modern  aircraft  are  raised  further,  and  with  the 
further  perfection  of  adaptive  control  theory  and 

the  development  of  computer  software  and  hardware  technologies,  the 
adaptive  scheme  could  expect  to  be  used  in  situations  where  the  airdata 
scheduling  scheme  does  not  work.  Accordingly,  the  advantages  of  the 
adaptive  control  scheme  such  as  saving  sensor  and  measurement  circuit¬ 
ry,  requiring  less  prior  knowledge,  saving  the  expenditure  and  time 
of  wind  test,  shortening  the  design  cycle,  providing  higher  flexibili¬ 
ty,  etc.  will  be  fully  utilized. 
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APPLICATION  OF  ROTOR  INDUCED  VELOCITY  FIELD 
OF  THE  VORTEX  THEORY  TO  THE  STUDY  OF 
STABILITY  AND  CONTROLLABILITY  OF 
HELICOPTER  WITH 
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( Air  Force  Research  Institute) 
Wang  Shicun 

(Nanjing  Aeronautical  Institute) 
Abstract 


l^th  the  aim  of  developing  a  new  rotor  system,  a  method  is  presen¬ 
ted  for  calculation  of  stability  and  controllability  of  helicopter  with 


hingeless  rotor.  The  effects  of  nonuniformity  of  rotor  induced  velocity 
distribution  and  blade  flapping  elasticity  are  taken  into  account  and  inve¬ 
stigated  in  detail. 

It  is  the  first  time  that  the  first  harmonic  term  of  rotor  induced  ve¬ 
locity  distribution  derived  from  the  generalized  vortex  theory  of  rotor 
is  employed  and  a  close  form  of  equation  for  the  induced  velocity  and  cir¬ 
culation  is  established.  Finally,  analytic  formulae  suitable  to  engineer¬ 
ing  application  are  given,  which  consider  the  effects  of  motion  parameter, 
control  parameter,  blade  flapping  elasticity  and  so  on.  Also,  the  effect 
of  the  second  flapping  mode  on  stability  and  controllability  of  helicopter 
with  hingeless  rotor  is  discussed  in  an  all-round  way.  The  calculations 
show  that  only  the  first  flapping  bending  mode  should  be  taken  into  acco¬ 
unt  in  stability  and  controllability  calculation  of  helicopter  with  hinge¬ 


less  rotor  in  the  range  of  ^<C0.3. 

Finally,  a  typical  helicopter  Bo-105  is  taken  as  a  numerical  example, 
and  it  is  found  that  tfie  results  are  in  good  agreement  with  its  flight-test 
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This  paper  presents  the  computation  method  of 
: 1  o  stability  and  controllability  of  a  helicopter 
hingeless  rotor.  The  effects  of  the  first 
harmonic  term  of  rotor  induced  velocity  distribu¬ 
tion  and  the  first  and  second  flapping  bending 
modes  of  a  rotor  blade  are  taken  into  considera¬ 
tion  and  investigated.  Especially,  the  applica¬ 
tion  of  rotor  induced  velocity  field  derived  from 
Kang  Shicun 's  generalized  vortex  theory  to  the 
study  of  stability  and  controllability  of  a  heli¬ 
copter  with  hingeless  rotor  is  described  in  detail. 
The  derivation  method  and  numerical  derivative 
method  for  Aly  and  are  employed.  Finally, 

the  Bo-105  helicopter  is  taken  as  a  numerical 

>mple,  and  the  computed  results  are  in  good  agree- 
"’•"Mt  with  its  flight -test  data. 


Nomenclature 


.  l  r. 

a0N'  a’N'  ®'>x  flapping  coefficients  corresponding  to  the  N1'  flapping 
mode  of  the  rotor  blade. 

( A) q ,  (A)1C»  (A)1S,  (A)°lc,  (A)°1S  Correction  coefficients  of  the 
monuniformity  of  rotor  induced  velocity  distribution. 

Mv„,  M„„,  mv„,  m„„  Rolling  and  pitching  rotor  hub  moments  and  corre- 
sponding  coefficients 
t  h 

S^(X^)  n  flapping  bending  vibration  model  function  of  a  rotor 
blade 

k  Number  of  rotor  blade 

T,  H,  S,  C  ,  C„,  C  Rotor  tension,  vertical  force,  side  force  and 

1  H  o 

corresponding  coefficients 
Rotor  composite  velocity  ( nondimensional  V^V^/riR) 

d.  1  Rotor  attack  angle  with  respect  to 

v^  Rotor  induced  velocity  ( nondimensional  v1=v^//2R) 

a  Rotory  rpm 

to)  ,  Rolling  and  pitching  angular  velocity  of  helicopter 

(  £>x=  cox//2 ,  £z=  *>z/o.  ) 

M,  My'  Mz  C^aracteristic  coefficients  of  flight  condition 

and  components  along  rotor  axis  system 

Frequency  ratio  corresponding  to  Nfc^  flapping  mode 

Generalized  coordinates  of  N^  flapping  mode 

Total  distance  angle 

Negative  torque  of  rotor  blade 
9,  Horizontal  periodic  distance  change  of  rotor 


K 

fr 


Vertical  periodic  distance  change  of  rotor 


Total  distance  angle  of  tail  rotor 


07TR 

P  Nonfactor  circulation  (f  =  r//2R2) 

I.  Preface 

The  rotor  induced  velocity  field  distribution  of  a  helicopter 
is  an  important  and  extremely  complicated  problem.  It  is  particularly 
important  to  a  helicopter  with  hingeless  rotor,  especially  the  first 
harmonic  term^1^.  In  the  past,  the  helicopter  balance,  stability 

and  controllability  computations  all  employed  the  assumption  of  uni- 

.  .  .  [1^-3  1 

form  distribution,  and  then  modified  methods  were  introduced  , 

yet  they  were  all  based  on  the  classical  momentum  theory.  These 
assumptions  are  too  simplified  such  that  not  only  the  effects  of 
lateral  disturbance  to  the  rotor  induced  velocity  field  are  not  con¬ 
sidered,  but  also  cannot  reflect  the  effects  of  various  structures, 
motion  and  control  parameter  of  a  helicopter.  The  development  of 
the  vortex  theory  is  becoming  more  and  more  mature,  and  using  it 
to  calculate  rotor  induced  velocity  field  distribution  is  more  accurate 
and  reasonable.  It  has  been  extensively  utilized  in  the  rotor  blade 
airloads  and  rotor  aerodynamics.  However,  due  to  the  complexity 
of  vortex  theory,  it  is  still  rarely  used  in  the  study  of  helicopter 
controllability.  This  paper  adopts  Wang  Shicun's  generalized  vortex 
theory  and  employs  the  close  form  of  equations  for  the  rotor  induced 
velocity  and  circulation  to  give  the  rotor  induced  velocity  field 
distribution  and  the  relational  equations  of  parameters  such  as  motion, 
structure  and  control  for  helicopter  with  hingeless  rotor  are  g;'ven. 
Moreover,  the  rotor  induced  velocity  field  of  the  vortex  theory  is 
introduced  into  all  the  hingeless  rotor  force,  moment  and  flapping 
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motion,  and  finally  analytic  formulae  suitable  to  engineering  applica¬ 
tion  are  obtained. 


The  effects  of  rotor  blade  elasticity  mode  number  to  the  features 
of  controllability  and  stability  for  a  helicopter  with  hingeless 
rotor  are  not  only  a  important  theoretical  problem,  but  also  a  practi 
cal  problem  in  engineering  application.  This  paper  has  made  detailed 
study  on  them. 


II.  Basic  Methods 


1.  Load  and  Velocity  at  the  Rotor  Blade  Element 


The  hingeless  rotor  blade  can  be  treated  as  a  suspended  beam. 
Also,  assume  the  blade  is  infinitely  rigid  in  twisting  and  chordvise 
bending  directions,  and  only  flexible  in  a  flapping  direction.  The 

blade  flapping  bending  differential  equation  is: 

aY  )  .  _  axY  aF 


(i) 


Jff 

mQ'rdr  is  the  centrifugal  tension,  dF/£r  is  the  airloads 


and  Co  r io 1 ' 


inertia  loads ; 


aE 

ar 


+  aF, 


ar  ar 

flow  region  and  compress ibi 1 i  tv  a  re  no* 


Here  the  effects  of  reverse 
considered. 


~r  — |-pO(QF)*(  +uI«r) 


where 


4>#  =  0,  -t-A^*+0iCO3»ti  +  9iSini)),  9,-«P7-0.7A9 


dr 


where 


The  velocity  at  rotor  blade  element: 

Sr**  *  +H,sini|>  +  H,cosiJ 

—  »o~  »,.C06\|>-5,,sini{i  +  (A,sint|>  +  Mtcostj>)  *  ^2) 

-  (^.costji  -  ^.siniji )  -  y/Q 

2 

;  y-Y/R,  y-  J  S„(*)U(<) 

AT-  1 

The  first  harmonic  distribution  of  the  rotor  induced  velocity 


is  considered. 


2.  Rotor  Induced  Velocity  Field 


According  to  the  blade  element  theory  and  Rukowski  formula, 
and  adopting  the  regular  assumption  that  U  *  Ux>  aw  Ql  ,  the 

circulation  relational  formula  is  obtained  as: 

f--|~a.b(«T<|»,+tt,)-f,+f,.cosiji  +  f1,sin^+ . 


Substituting  u_,  u  and  <£  -  and  after 

i  p 

m 

f.— i-a.b[9,*+A<P** 

+ ) 
f1<-^-o.b[01*+M,+li.A<P*-51,+c3^f+2sK6lN 
els#  1 

-*•***-&  J 

f  „  » -y-o.b  [etx + ^.0, + -  V,,  +  cb.X  -  SSjvfljA 

+^2a*K4r] 

Ac'ord-.g  to  the  stationary  vortex  theory  f 


[4] 

'lucun  ,  and 


only  limited  to  the  first  harmonic  term,  the  rela^i  al  formula  between 

i  r.  :®d  velocity  and  circulation  is: 


5 _ fe^f-Jcosa^/  _2_+x+_l_.xt\  +_2snmJ_f<'| 

**  4KV]  l  1  +sii>a,  \  3  2/1  +sina,  J 

— * — f — 2COS.a‘ — y,(  3  -3x+**)f,H — .  2  -  -f„  1 

1  4«Vi  l  1  +sina,  YlK  1  +stna,  J 


Solving  (3)  and  (4)  simulataneously  and  eventually  the  induced 

velocities  are  obtained:  B.-5.+5..«**+5.*«i»4  + . 

V, «  (  A  )  ,  [e^+  H.0.+ -y-H.  0, + M:n)  -37-] 

5,.-  (y4)..[e,*+»1.0,+  ^  A?*+Spr+2(*A*-y.a,N-^Jj 

+u)  iJ(— f-+  * +4-*1)[®^+A*^-4'+4-^ei+^Mi 
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where  .t  definitions  of  (A,'„,  (A)^,  (A)°c,  (A)^s  an<^  iS  are 
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3.  Flapping  Motion 


Using  the  rotor  blade  vibration  model  equation  and  the  model's 
normal  feature,  then  from  the  blade  flapping  bending  equation  the 
rotor  flapping  motion  equation  can  be  obtained: 


gUj) 


+  «W4) 


IV  =*  1  ,  2 


(6) 


Generally  speaking,  it  is  difficult  to  accurately  solve  the 
vibration  model  equation.  This  paper  adopts  the  "asymptotic  method" 
to  obtain  the  solution.  It  is  already  sufficient  to  select  the  first 


harmonic  term. 


C.v(  ^  )  =a9l4  —  a;AcosiJi—  6,.vSiTiij> 


then  the  above  equation 


can  be  changed  to: 


[ 


A*J0.v-(X£-  l  HnCOS\|>  — Ol—  1  )6,/.sintJ> 


l  mSXx)dx 


In  essence,  Equation  (6)  and  (7)  describe  the  forced  effects 
of  each  order  vibration  model  S^(x)  under  external  loads.  As  a  cen¬ 
trally  hinged  type  rigid  blade,  ^^  =  1,  (x)-x,- and  Equation  (7) 

_. ,  rR  dF  j 

is  simplified  to  the  hinged  type  rotor  flapping  equation  o“5Trdr’ 

with  which  we  are  familiar.  Obviously,  the  hinged  type  rotor  blade 
is  a  special  case  of  hingeless  rotor.  According  to  the  aforementioned 
external  loads  and  induced  velocity  distribution,  and  eliminating 
harmonic  terms  above  the  second  order  and  those  cross-coupling  terms 
between  each  order  mode,  two  groups  of  linear  algebraic  equations 
can  be  eventually  obtained: 


(81 


CCHW-C'DJ 

correspond i ng  first  order  node:  Cs 0  =■  (a(,aM6,|]r 


corresponding  second  order  mode: 


C>»iD = 


When  flight  conditions  are  given,  [A],  [  B ]  ,  [C]  and  [D]  are 

all  determined.  The  [  and  [  £  2 1  can  be  solved  and  the  correspond¬ 

ing  derivatives  can  be  obtained. 


4.  Hingeless  Rotor  Force  and  Rotor  Hub  Moments 

The  rotor  hub  moments  are  particularly  important  to  the  hingeless 
rotor.  Take  the  integral  of  moment  generated  by  the  aerodynamic 


force  on  the  blade,  cetrifugal  force  and  the  flapping  and 
Coriolis  inertia  forces  with  respect  to  the  center  of  rotor: 


Af„(  0  .  t  )*  j  ^(my+mQ:K — *Jr~)ri/r 


If  only  the  first  order  flapping  coefficient  is 
rotor  hub  moments  generated  by  the  centrifugal  force 
inertia  force  cancel  each  other  exactly.  Therefore, 

M - Trf'.' 

M-  —  fr  .T,"  I » ’XT  ' 


selected,  the 
and  the  flapping 


Take  rotor  blade  elasticity  first  and  second  flapping  modes 
and  first  harmonic  term  of  induced  velocity  into  account,  then  the 
following  can  be  obtained  after  integration  and  nondimensionalization : 

m,„  »  m*«0  ■+■  +  Am.* 

Similar  to  a  hinged  rotor,  the  elementary  forces  of  a  hingeless  rotor 
are:  dT^dL**—-  Pa»6(«|<j)#+uTuf)c/r 

dH  S?  (  d  D — rfZ-3,  )sin  i|i  —  dL  cos^ 

dS^  —  (dD  —  dL$*)cos'b—dL-^^-sm'i} 


After  the  same  man i f upla t ion ,  the  following  can  be  obtained: 

Cj“Cr,  +  ACr»  ChjH" AChj  c,«c„+ac, 

CTQ,  r  :s0,  mxH0  and  mzHQ  are  the  basic  terms;  A  C^,  A  C^,  ACg, 
A  m  and  A  m^^  are  the  increments  derived  from  induced  velocity 

harmonic  terms  in  which  this  paper  is  very  interested'. 


i  J  0  6  J  0 

A CH  =  -y-  0,  J  A0>  J  *  **5,,  djf  +  -T  05“<f*  +  ,1,J  ol 

H — 8.(4  H.j  o«V  </*-  J#l'ot'i/^ 

+ 4-^ . .[«, ,  J  \  dx+ alt  J  J  + -r  M*-  J  0  D“ 'sr 

+  6>‘I  o5*'"^*] — r^W  ldit^dx+  fluf  o5,‘^H*] 
— r*'[buS \i^Js+buS  0  **&*] 
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According  to  the  induced  velocity  distribution  of  the  vortex 
theory,  the  hingeless  rotor  force  and  rotor  hub  moment  which  take 
the  induced  velocity  harmonic  components  and  blade  elasticity  into 


(  14) 


account  can  bn  eventually  obtained: 

— —  C F 1  -  Cfl.3  < E  >  +  C A/?)  <E  > 

a_o 


where 


CF}  — CCr  CH  C,  m,H 

{£}«<0,  A*P  4,  9,  9,  »,  4,  1  )T 

[Rq]  and  [AR]  are  functions  of  and  (A)Q,  (A)1C»  (A)°c,  (A)lg, 

(A)°s  is  the  flight  conditions  and  control  parameter  matrix. 

Similarly,  when  X|*  1 ,  St(x)~  xt  Xi**0,  —  0  ,  the  ordinary  hinged 

rotor  forces  and  rotor  hub  moments  can  be  obtained. 

5.  Disturbance  Motion  Equations  and  Control  Derivatives 

The  linear  small  disturbance  motion  differential  equation  groups 
are  divided  into  longitudinal  and  lateral  coupling  and  noncoupling 
for  consideration.  The  disturbance  motion  equaiton  is: 

an<H>-cin{G>  (15) 

where  [W]  is  the  stability  derivative  matrix;  [U]  is  the  control 
derivative  matrix;  [hJ  is  the  motion  parameter  matrix;  and  {g]  is 
the  control  parameter  matrix. 

The  coefficients  of  the  equation,  i.e.,  all  the  stability  deriva¬ 
tives  in  matrix  [ W ]  and  the  control  derivatives  in  matrix  [U],  are 
obtained  by  computer  using  the  triple-point  formula  numerical  deriva¬ 
tive  method  with  step  increment  H=0.005.  The  computations  have  shown 
that,  for  stability  and  controllability,  the  numerical  derivative 
method  is  feasible  and  conducive  to  the  rotor  aerodynamic  and  dynamic 
model's  being  more  truthful  and  perfect. 


j  . 


III.  Examples 


In  order  to  explain  the  method  of  this  paper,  a  Bo-105  helicopter 
with  hingeless  rotor  is  selected  as  a  numerical  example  and  the  fuse¬ 
lage  data  is  substituted  by  test  results  of  similar  fuselage  (longi¬ 
tudinal  data  are  taken  from  Reference  [5],  lateral  data  are  taken 
from  H-5).  See  Reference  [6]  for  original  data  and  the  balance  data 
are  taken  from  Reference  [7],  The  "asymptotic  method"  is  adopted 
for  the  blade  frequency  vibration  model  to  obtain  ^  =  1.124,  \2  =  2.86. 
Four  kinds  of  flight  conditions  at  0.1,  0.2,  0.3  for  hovering 
and  cruise  are  calculated  under  three  situations  of  considering  both 
the  induced  velocity  harmonic  component  and  second  flapping  mode 
and  with  either  one  eliminated.  In  order  for  convenient  comparison, 
Fig.  1.  shows  a*1  =  a^-99;  b*  ]_  =bi  ]_  ~e  •  It  can  be  seen  that  the  induced 
velocity  harmonic  term  primarily  affects  the  lateral  flapping  coeffi¬ 
cients.  There  is  a  peak  at  /A  ~ 0 . 1  and  this  phenomenon  has  been 
observed  in  many  experiments.  This  result  matches  completely  with 
those  of  Reference  [3]  and  [4].  It  can  be  seen  from  Fig.  2  that, 
within  the  range  of  < 0.3,  the  effects  of  the  second  flapping  do 
not  excee  .  10%  of  the  total  value.  When  ft- >  0.3,  the  effects  of 
the  second  mode  must  be  considered.  The  analysis  of  this  paper  is 
basically  consistent  with  the  conclusions  of  Reference  [1]  and  [8]. 

The  results  have  shown  that  the  vortex  theory  induced  velocity  dis¬ 
tribution  introduced  by  this  paper  is  reasonable  and  correct. 


(7)  extended  line;  (8)  first  mode;  (9)  total  moment;  (10)  second 
mode;  (11)  total  moment. 


The  stability  and  controllability  features  under  various  condi¬ 
tions  are  computed.  Figure  3  shows  the  comparison  of  computed  result 
with  the  theoretical  values  and  flight-test  data.  Also,  the  primary 
stability  and  control  derivatives  for  typical  flight  conditions  are 
selected  to  compare  with  the  theoretical  values  and  flight  results 
given  b‘  Reference  [9],  and  they  are  in  better  agreement. 


Fig.  -.  '  mparison  of  characteristic  root  resul;.-. 

Key:  (1;  Reference  [9];  (2)  flight-test  data  of  eohod;  (3)  flight 

test  data  of  ML  method;  (4)  flight-test  data  of  ML  method  (nonlinear 
model);  (5)  quasi-steady  theory  computed  values;  (6)  computed  results 
of  this  paper;  (7)  third  class. 


Table  1 . 

Key:  (1)  derivative;  (2)  This  paper;  (3)  unit;  (4)  Computed  results 

of  this  paper  1750  kg,  43.6  m/sec;  (5)  All  taken  into  account;  (6) 
data  of  Reference  [9]  2100  kg,  38.1  m/sec;  (7)  quasi-steady  theory; 
(8)  computed  values;  (9)  flight-test  data;  (10)  minimum  second  multi¬ 
plication  method;  (11)  maximum  approximation  method. 


IV.  Conclusions 


This  paper  establishes  for  the  first  time  the  stability  and 
controllability  computations  of  a  helicopter  with  hingeless  rotor 
while  taking  the  rotor  induced  velocity  first  harmonic  distribution 
of  vortex  theory  into  account.  It  presents  analytic  formulae  suitable 
to  engineering  applications.  Within  the  range  of  /“<  0.3,  the  calcula¬ 
tion  only  needs  to  consider  the  first  flapping  mode. 


The  example  computations  of  this  paper  were  primarily  completed 
by  Sergeant  Engineer  Chen  of  Department  602;  engineers  Miar  Qiuan, 

Cao  Gebing  and  Ho  Rubiang  also  provided  their  assistance.  During 
the  research  process,  we  received  guidance  and  assistance  from 
Professor  H.  C.  Curtis  of  Princeton  University  and  Associate  Professor 
Xu  Xingyu  and  Instructor  Chen  Qiuming  of  Nanjing  Aeronautical  Institute. 
We  hereby  express  our  thanks  to  all  of  them. 
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A  SECOND-ORDER  THEORY  FOR  THREE- 
DIMENSIONAL  UNSTEADY  FLOWS  AND 
ITS  APPLICATION  TO  THIN  WINGS 
WITH  ANGLE  OF  ATTACK  AT 
SUPERSONIC  AND  HYPE-  ‘ 

RSONIC  SPEEDS 

He  Longde,  Han  Yanliang,  Quirt  Fuxing 
( Institute  of  Mechanics,  Academia  Sinica ) 

Abstract 

^^ccording  to  the  assumption  of  high  Mach  numbers  the  secondorder 
equation  for  unsteady  flow  reduces  to  a  form  analogous  to  that  for  steady 
flow.  As  an  approximation)  a  particular  solution  is  derived  from  the  assu¬ 
mption  of  local  two-dimensionality  instead  of  the  approximate  three-di¬ 
mensional  particular  solution. 

In  treatment  of  the  unsteady  problem  for  delta  wings  with  low  aspect 
ratio  and  supersonic  leading  edges,  the  non-linear  pressure  relation  and 
the  modified  reference  flow  parameters  are  adopted  to  take  the  effect  of 
mean  angles  of  attack  into  account.  Analytical  expressions  are  obtained  for 
the  first  order  expansion  of  the  unsteady  velocity  potential.  They  condu¬ 
ce  to  simplifying  the  calculation  of  stability  derivatives  to  a  great  ex- 


Stability  derivatives  for  flat  plate,  double-wedge  or  biconvex  delta 
wings  oscillating  at  low  frequencies  are  given  respectively.  Since  the  cor¬ 
responding  experimental  data  are  unavailable,  the  obtained  results  are 
merely  compared  with  thossjb!  Hui's  flat  pIateC4J,  and  both  results  are  in 
good  agreement.  ^ '  ' 


A  SECOND-ORDER  THEORY  FOR  THREE-DIMENSIONAL  UNSTEADY  FLOWS  AND  ITS 
APPLICATION  TO  THIN  WINGS  WITH  ANGLE  OF  ATTACK  AT  SUPERSONIC  AND 
HYPERSONIC  SPEEDS 

He  Longde,  Han  Yanliang  and  Qian  Fuxing 
Institute  of  Mechanics,  Academia  Sinica 

I.  Preface 

The  dynamic  derivatives  are  indispensable  parameters  in  the 
dynamic  quality  analysis  of  aircraft  and  the  design  of  its  control 
systems;  they  directly  affect  the  dynamic  characteristics  of  aircraft. 
Therefore,  since  the  70'-s,  there  have  been  many  studies  made  either 
in  wind  tunnel  tests,  flight  tests'  or  analyses',  computation,  etc., 
and  more  progress  has  been  made^1^. 

The  computation  of  dynamic  derivatives  are  essentially  unsteady 
aerodynamic  computation  problems.  The  existing  linear  theories  cannot 
satisfy  the  requirements,  and  the  demand  for  taking  nonlinear  effects 
into  account  in  the  dynamic  derivative  computation  has  grown  more 
urgent  by  the  day.  The  second-order  theory  is  one  kind  of  method 
to  compute  nonlinear  effects.  It  is  established  by  second  iteration 
or  the  second-order  expansion  of  the  object  surface  deviation  angle, 
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and  the  primary  nonlinear  effects  are  computed  from  the  second-order 


equation.  Van  DykeL  J  and  Qian  Fuxing  have  worked  in  this  area. 


Whereas  when  subtraction  is  performed  between  the  pressure  coeffi¬ 


cients,  which  have  undergone  the  second-order  theory  series  expansion, 


of  the  upper  and  lower  wing  surfaces,  the  unsteady  effects  of  the 


angle  of  attack  are  cancelled  out;  therefore.  Reference  [3]  only 


computed  the  pitching  damping  features  of  the  compression  surface 

-  .  .  r  4 1 

of  a  three-dimensional  thin  wing.  In  recent  years,  Hui  et  al. 


have  proposed  a  unified  supersonic  and  hypersonic  theory  which  can 


give  better  results  for  low  frequency  flat  plate  wing  with  high  Mach 


number  and  angle  of  attack,  but  it  is  not  suitable  to  wings  with 


a  small  aspect  ratio  under  low  supersonic  conditions,  and  it  cannot 


handle  thickness  effects. 


This  paper,  on  the  basis  of  a  supersonic  unsteady  three-dimension¬ 


al  linear  solution,  utilizes  the  approximate  local  two-dimensional 


results  of  the  second-order  theory  at  a  high  Mach  number  to  take 


the  mean  angle  of  attack  and  nonlinear  effects  of  thickness  into 


account  and  specifically  solves  the  unsteady  problem  of  supersonic 


delta  wings  with  angle  of  attack  at  leading  edges.  In  order  to  compute 


the  nonlinear  effects  of  the  angle  of  attack,  this  paper  conducts 


computations  using  the  nonlinear  isothermal  pressure  relation  and 


the  reference  incident  flow  parameters,  and  compares  with  the  results 


of  Reference  [4].  They  are  in  good  agreement. 


The  basic  equation  of  the  unsteady  second-order  theory  and  its 


solution  are  omitted  here.  See  Reference  [3]  for  details.  It  needs 


to  be  explained  that,  when  solving  the  unsteady  second-order  equation, 


our  method  adopts  the  assumption  of  a  high  Mach  number  approximation 
and  local  two-dimensionality  of  the  flow  on  a  thin  wing  surface. 

II.  Nonlinear  Isothermal  Pressure  Relation  (Method  1) 

A  three-dimensional  thin  wing  with  a  small  aspect  ratio  in  a 
uniform  supersonic  flow  is  subject  to  small  amplitude  simple  harmonic 
pitching  oscillation  near  the  mean  angle  of  attack  If  second- 

order  is  considered,  the  flow  is  isothermal.  Rectangular  coordinates 
oxyz  are  selected  with  the  origin  at  the  tip  of  the  wing's  leading 
edge,  x-axis  along  the  flow  direction,  y-axis  pointing  to  the  right 
when  facing  the  flow  and  z-axis  being  vertical  with  up  as  positive, 
t  represents  time.  Assume  <#>(x,  y,  z,  t)  is  the  disturbance  velocity 
potential,  then  the  formula  for  nonlinear  pressure  coefficient  is: 

]  V"  1  “  l}  (2.1) 

where  ,  Ma &  >  f  and  H  are  free  flow  velocity,  Mach  number,  adia¬ 
batic  index  and  height  of  unsteady  object  surface,  respectively; 
the  subscript  respresents  a  deviation  derivative  such  as  <P=  J<P/J  x. 
(p  can  be  represented,  using  the  sum  of  steady,  unsteady  first  and 
second  disturbance  velocity  potentials,  as^^ 

<&  =■<»>.(*,  y,  *  )+*,(*,  y, *  *,  y,  z) 

+*«<*,  y,  on  1 

By  expanding  Equation  (2.1)  and  omitting  a  small  quantity  of  higher 
order,  the  pressure  coefficient  formula  of  the  second-order  theory 
is  obtai ned 

C,u -  -  ( «. +  'K «-»•«) 

+  *  +«#*)-«*'(  *-b  >*•«) 


where  i  =»/-!,  jc  =  <D/r.,  o—&Mol/(Mai-  1  )»  CO  and  9n  are  the  Angular 


frequency  and  angular  amplitude  of  the  simple  harmonic  oscillation; 


h  and  b  are  the  height  of  steady  object  surface  and  x  coordinate 


of  the  vibration  center.  tfo'  4*0  ars  t^ie  loca-!-  two-dimensional 

solutions,  ip*  is  the  unsteady  three-dimensional  linear  solution 


without  the  high  Mach  number  approximation  to  take  the  three-dimen¬ 


sional  effects  into  account  and  raise  the  computation  accuracy  under 


low  supersonic  conditions. 


When  computing  the  dynamic  derivatives,  the  pressure  difference 


Cpu~Cpi  between  the  upper  and  lower  surface  will  be  integrated 


along  the  . object  surface.  If  the  upper  and  lower  wing  surfaces 


are  symmetrical,  the  difference  between  the  solutions  of  upper  and 


lower  surfaces  is  only  the  sign  change  of  (X  n  and  9n.  It  can  be 


seen  from  Equation  (2.3)  that,  in  the  second-order  theory,  there 


is  no  9q  cX  q  effect  included  in  Cnil-Cn1  .  Therefore,  Equation  (2.3) 


pu  pi' 


cannot  be  used  to  obtain  the  unsteady  effect  9^  (X  q  of  the  angle  of 


attack  using  the  velocity  potential  of  the  second-order  theory;  the 


nonlinear  Equation  (2.1)  of  the  pressure  coefficient  should  be  used. 


First,  convert  Equation  (2.1)  to  obtain  solution  at  z=0,  then  solve 


for  derivative  with  respect  to  9,  and  take  limit  9Q-*  0  and  omit  a 


small  quantity  of  higher  orders,  thus 


“5<f C'  I*-*-  c  1  -(  Y  -  1  Y  - 1  c; 


'2.-1) 


i  to  t, 

where  9  =  9qC  ;  the  representation  of  C*  is  similar  to  Equation 


(2.3),  i.e.,  Equation  (2.3)  divided  by  9  and  then  adding  -2e 


-  ikx2 


‘  A'VvVi’ 


The  unsteady  pitching  moment  coefficient  on  the  compression 


surface  or  expansion  surface  is 


lim  = -J- J  J  (  *  - 

9,~o  *8  SciJ 


b  )C  1  -(  Y  -l)A/ai<P„)  v  ~  i  C;dxdy 


(2.3; 


1  CC«#  +  ihCm  j)  .■ 


where  C  £  and  CmQ  are  the  pitching  damping  derivative  and  static 


'mQ 


derivative;  9=  (  dQ/dt )  •  (  Cq/V^q  )  =  ik9Qe 1  ,  k  =COc0/V00  is  the  folding 

frequency;  c^  is  the  chord  length  at  wing  root;  S  is  the  surface 
area  of  wing.  For  low  frequency,  (p  is  cal.cul.ated  to  k^  term  using 
the  frequency  series  expansion  method.  Here  "one  surface"  means 
either  the  compression  surface  or  the  expansion  surface,  and  the 
total  stability  derivative  is  the  algebraic  sum  of  upper  and  lower 
surfaces.  The  sign  of  (X  q  must  be  changed  on  the  upper  and  lower 
surfaces  of  the  wing.  In  actual  computation,  the  wing  crosssection 
considered  is  near-jet  similar  with  the  object  surface  equations 
as : 

double-wedge:  2(5//)-  e  Cc.-|c,-2*/ / 1}/ 2 

double-arc  :  Z(5/  /  )"2  «(*//)•(  1 

where  ;c=x-ytgXs;  f  =  i -y  tgXs/Cg ;  Xg  is  the  swept  angle  of  the  wing's 
leading  edges;  flat  plate  corresponds  to  the  situation  of  thickness 
ratio  &  =0 . 


III.  Reference  Incident  Flow  Parameter  (Method  2) 


Based  on  the  second-order  theory  of  isothermal  flow,  the  change 
in  enthalpy  after  the  flow  crosses  the  shocks  wave  is  completely 
ignored.  Thus,  under  the  conditions  of  a  high  Mach  number  and  larger 
angle  of  attack,  the  error  in  the  second-order  theory  increases 
as  the  shock  wave  strengthens.  The  reference  incident  flow  parameter 
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method  is  exactly  a  method  proposed  to  solve  the  problem  of  computing 
steady  effects  of  the  angle  of  attack  and  thickness  effects. 


When  a  two-dimensional  steady  supersonic  flow  flows  over 


a  flat  plate  at  angle  of  attack  # g ,  the  undisturbed  free  flow  param¬ 


eters  of  the  streamlines  of  a  three-dimensional  thin  wing  are  replaced 
by  flow  parameters  of  the  object  surface  behind  the  shock  wave/  ■ 
then  the  upper  and  lower  surfaces  of  the  wing  are  computed  respective¬ 
ly  using  the  second-order  theory.  The  reference  incident  flow  param¬ 
eters  can  be  completely  determined  by  the  shock  wave-expansion  wave 
method  and  is  denoted  by  a  subscript  1.  It  should  be  pointed  out 
that  in  Method  2  the  mean  angle  of  attack  C does  not  appear  directly 
and  its  effects  are  included  in  the  reference  incident  flow  param¬ 
eters.  The  formula  for  pressure  coefficient  is 


P-P-  P-Pi  , 


1 


-P.Vl 


1 


-P.F  i 


-PmVl 


(3.1) 


2  -  2  '  2 

The  second  term  to  the  right  of  the  equality  sign  is  the  steady  pres¬ 


sure.  It  is  independent  of  the  stability  derivatives  Cm£  and  Cmg 


and  thus  is  not  considered. 
Let  — 


p  — a  _  p y\  .  p  -Pi 

1  _ o  y »  P-^i  1 

2 


-~pmvi  r-T "  -—pyi 


G(Mam,  ac)CJ 


(3.2) 


G(Ma„  at)-p^>/p.FI 


(3.3) 


where  p ,p  and  V  are  pressure,  density  and  velocity,  respectively. 
Using  Ma*o  and  $q,  the  G  of  the  compression  surface  and  expansion 


surface  can  be  obtained.  The  form  of  (V  is  identical  to  Equation 


(2.3),  except  that  the  undisturbed  flow  parameters  are  replaced 
by  the  reference  incident  flow  coefficients.  Take  the  integral 


of  Cp  on  the  wing,  then  the  unsteady  aerodynamic  forces  and  the 
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desired  stability  derivatives  are  obtained. 

IV.  Results  and  Discussions 

It  can  be  observed  from  Figs.  1  and  2  that,  for  a  flat  plate 
delta  wing,  the  results  of  this  method  and  the  theoretical  results 
of  Reference  [4]  are  in  good  agreement  both  in  the  trend  and  order 
of  magnitude.  '  "  • 

For  a  delta  wing  with  double-wedge  crosssection,  the  results 
of  the  two  methods  in  this  paper  have  shown  that,  when  Ma*x>  d  q<  1 . 3 , 
the  consistency  of  the  two  methods  are  better.  Under  high  Mach 
number  condition.  Method  2  is  more  reliable.  However,  at  Ma*x><2.5 
and  larger  cL  Q/  the  reference  incident  flow  Mach  number  Ma^^  on  the 
compression  surface  is  smaller  and  the  nearsonic  or  subsonic  leading 
edges  will  occur.  At  such  time,  Method  2  becomes  useless  and  Method 
1  can  be  used.  The  two  methods  are  complementary  to  each  other. 

They  can,  within  the  range  of  Mach  number  2  8  and  folding  frequency 

0 ~ 1 . 0 ,  more  accurately  estimate  the  thickness  effects,  especially 
the  unsteady  effects  of  angle  of  attack. 

This  method  is  both  analytic  and  semi-analytic.  It  is  convenient 
to  use  and  its  operation  fast.  It  can  combine  the  unsteady  second- 
order  theory  results  of  a  single  wing  and  the  theoretical  interference 
factors  of  a  long  thin  object  to  conduct  unsteady  aerodynamic  computa¬ 
tions  for  the  assembled  body  of  wings,  fuselage  and  tail  wings. 
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Fici.  ..  Variation  of  with  Mach  number  Ma*#- 

Kev  :  delta  wing;  (2)  flat  plate;  (3)  Method  _  this  paper; 

(4)  Reference  [4];  (5)  double-wedge;  (6)  Method  2  of  this  paper; 
(7)  Method  1  of  this  paper. 


Fig.  4.  Variation  of  Cmg ,  with  ••••.'  of  oscillation  b/Cg 

Key;  (1)  delta  wing;  (2)  flat  plate;  (3)  Method  2  of  this  paper; 

(4)  Reference  [4];  (5)  double-arc;  (6)  Method  2  of  this  paper. 

At  a  low  supersonic  speed  such  as  Ma^o  =2,  the  accuracy  of  the 
second-order  solution  is  not  high.  However,  since  one  uses  the 


which  does  not  have  as 


unsteady  three-dimensional  solution  *  , 
high  Mach  number  approximation,  the  overall  computation  accuracy 
can  be  assured,  and  the  linear  results  are  improved. 
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